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We re-discuss the evolutionary state of upper main sequence magnetic stars using a sample of Ap and Bp stars with 
accurate Hipparcos parallaxes and definitely determined longitudinal magnetic fields. We confirm our previous results 
obtained from the study of Ap and Bp stars with accurate measurements of the mean magnetic field modulus and mean 
quadratic magnetic fields that magnetic stars of mass M < 3 Mq are concentrated towards the centre of the main-sequence 
band. In contrast, stars with masses M > 3 M© seem to be concentrated closer to the ZAMS. The study of a few known 
members of nearby open clusters with accurate Hipparcos parallaxes confirms these conclusions. Stronger magnetic fields 
tend to be found in hotter, younger and more massive stars, as well as in stars with shorter rotation periods. The longest 
rotation periods are found only in stars which spent already more than 40% of their main sequence life, in the mass 
domain between 1.8 and 3Mq and with logg values ranging from 3.80 to 4.13. No evidence is found for any loss of 
angular momentum during the main-sequence life. The magnetic flux remains constant over the stellar life time on the 
main sequence. An excess of stars with large obliquities f3 is detected in both higher and lower mass stars. It is quite 
possible that the angle f3 becomes close to 0° in slower rotating stars of mass M > 3 M© too, analog to the behaviour 
of angles /3 in slowly rotating stars of M < 3 M©. The obliquity angle distribution as inferred from the distribution of 
r-values appears random at the time magnetic stars become observable on the H-R diagram. After quite a short time spent 
on the main sequence, the obliquity angle (3 tends to reach values close to either 90° or 0° for M < 3 M© . The evolution 
of the obliquity angle j3 seems to be somewhat different for low and high mass stars. While we find a strong hint for an 
increase of /3 with the elapsed time on the main sequence for stars with M > 3M©, no similar trend is found for stars 
with M < 3 M©. However, the predominance of high values of (3 at advanced ages in these stars is notable. As the physics 
governing the processes taking place in magnetised atmospheres remains poorly understood, magnetic field properties 
have to be considered in the framework of dynamo or fossil field theories. 
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!> . 1 Introduction studies of the evolution of some characteristics of magnetic 

Ap and Bp stars were based on smaller stellar samples (e.g., 
j^j ■ After presenting in the first part of this paper the results of a North 1985, Hubrig, North & Mathys 2000, Hubrig, North 
\ magnetic field study of more than one hundred chemically & Szeifert 2005), or on stellar samples which included CP 
peculiar (CP) A and B stars with previously unknown or stars with not well-defined magnetic field strengths or rota- 
only poorly determined magnetic fields using FORS 1 at the tion periods (e.g., Rudiger & Scholz 1988, Pohnl, Paunzen 
Very Large Telescope (Hubrig et al. 2006), we examine in & Maitzen 2005, Kochukhov & Bagnulo 2006). 
this second part of our paper the evolution of the magnetic 
field across the upper main sequence (MS). The emphasis is 
put on the study of the evolutionary aspect of the magnetic 
field geometry in stars with accurately determined positions 
in the H-R diagram and well-studied periodic magnetic vari- 



The determination of rotation periods of Ap and Bp stars 
is in some cases not an easy task, as some stars are ex- 
tremely slowly rotating with periods up to 70 years and 
longer. The fact that magnetic stars are slow rotators sug- 



„ , , c , . - c c i j gests that some mechanism of magnetic braking has been 

attons. Knowledge of the evolution of the magnetic field ° ° ° 

,. " . „ ,. . ... working in a certain phase of the stellar evolution. However, 

and its geometry, especially of the distribution of the obhq- ° r 

, a , . . .. . • -4.1. . the comparison of rotation periods of Ap stars of different 

uity angle a (orientation of the magnetic axis with respect * T 

. . .. ... .. , . j . j , . , ages did not present any evidence supporting the hypothesis 

to the rotation axis) is essential to understand the physical ° cc a . . ,. , . 



processes taking place in these stars and the origin of their 
magnetic fields (fossil versus contemporary dynamo). Early 



that Ap stars suffer considerable magnetic braking during 
their MS life (e.g., North 1985, North 1998, Hubrig, North 
& Mathys 2000). Also the search for progenitors of slowly 
Corresponding author: shubrig@eso.org rotating magnetic Ap stars among normal A stars showed 

no statistically significant difference between the v sin i dis- 
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tributions of young and old A-type stars (Hubrig, North & 
Medici 2000). 

The periodic magnetic variations of Ap and Bp stars 
are generally described by the oblique rotator model (Stibbs 
1950) with a magnetic field having an axis of symmetry at 
an angle to the rotation axis, usually called (3. If the mag- 
netic field is excited by a dynamo mechanism, the magnetic 
field can be either symmetric or antisymmetric in regard to 
the equatorial plane (Krause & Oetken 1976). In case the 
stars have acquired their magnetic field at the time of their 
formation or early in their evolution, the currently observed 
magnetic field is called a fossil field. Randomness in the 
(3- distribution is usually regarded as an argument in favour 
of the fossil theory, as the star-to-star variations in obliq- 
uity of magnetic field axes can plausibly be interpreted as 
reflecting differences in the intrinsic magnetic conditions 
at different formation sites. Clearly, the decision between 
the fossil and dynamo explanations for the magnetic field 
origin can be made only through careful consolidation of 
magneto-hydrodynamic theory and observations. Our find- 
ing that no strong magnetic fields are found in A stars close 
to the ZAMS with masses M < 3 M©, supports those theo- 
ries that locate a dynamo in the convective core since there 
should be enough time in the star's evolution to transport 
the core dynamo field to the stellar surface. As hotter, more 
massive magnetic Bp stars are mainly found to be located 
on the ZAMS or very close to it (Hubrig, North & Szeifert 
2005) it is quite possible that Bp stars acquired their mag- 
netic field already at the time of their formation. The previ- 
ous weakness of the fossil-field theory has been the absence 
of magnetic field configurations stable enough to survive 
in a star over its lifetime. Braithwaite & Spruit (2004) re- 
cently carried out numerical simulations to show that stable 
magnetic field configurations can develop through evolution 
from arbitrary, unstable initial fields. The possibility of a 
magnetic field resulting from a magneto-rotational instabil- 
ity in the radiative envelope has been discussed by Riidiger, 
Arlt & Hollerbach (2001) and Arlt (2004). Over the years, 
a lot of work has been done to show that differential ro- 
tation, advection by internal circulation, stellar wind, mag- 
netic torques, and rotation about a non-principal axis can 
all produce changes in the (3 angle of an oblique rotator. 
Although the time scales over which these (and other) pro- 
cesses cause (3 to change remain fairly uncertain, it could 
be possible that the observed magnetic field distributions 
reflect not only how the magnetic field was produced but its 
subsequent evolution in response to the operation of various 
mechanisms. 

2 Basic data 

To study the evolution of the magnetic field distribution 
in Ap and Bp stars we decided to compile Hipparcos par- 
allaxes and photometric data exclusively for magnetic Ap 
and Bp stars with well studied periodic variations of their 
mean longitudinal magnetic fields. The mean longitudinal 



magnetic field is an average over the visible stellar hemi- 
sphere of the component of the magnetic vector along the 
line of sight, and is derived from measurements of wave- 
length shifts of spectral lines between right and left cir- 
cular polarisation. The individual published measurements 
of the mean longitudinal magnetic field (B{) have recently 
been put together in the catalog of stellar magnetic rota- 
tional phase curves (Bi) (<fi) by Bychkov, Bychkova & Madej| 
(2005). The phase <f) is determined by 

4> = 2*(^±), (1) 

with Ti the time of measurement, T the zero epoch, i.e. the 
time corresponding to the phase <f> = 0, and P the rotation 
period. Out of 127 Ap and Bp stars presented in this catalog, 
only for 89 stars with well defined magnetic curves the par- 
allax error is less than 20% and Stromgren and/or Geneva 
photometry is available. One more star, HD 81009, with an 
accurate Hipparcos parallax and a well-defined magnetic ro- 
tational phase curve studied by Landstreet & Mathys (2000) 
has also been included in our sample. 

Our results of the study of the evolutionary state of up- 
per main sequence magnetic stars using a smaller sample 
of Ap and Bp stars with accurate Hipparcos parallaxes and 
definitely determined longitudinal magnetic fields indicated 
a notable difference between the distributions of Ap and Bp 
stars in the H-R diagram. In contrast to magnetic stars of 
mass M < 3M which have mostly been found around 
the centre of the main-sequence band, stars with masses 
M > 3 M seemed to be concentrated closer to the ZAMS, 
and the stronger magnetic fields tend to be found in hotter, 
younger (in terms of the elapsed fraction of main-sequence 
life) and more massive stars (Hubrig, North & Szeifert 2005).| 
In view of these apparent differences in the evolutionary 
state between magnetic stars in different mass domains, we 
decided to group our sample in two sub-samples with one 
containing lower mass stars with M < 3 M Q (57 stars in 
all) and another one with M > 3 M (33 stars). The ef- 
fective temperatures of the stars of both sub-samples have 
been determined from photometric data in the Geneva sys- 
tem through the calibration of Kiinzli et al. (1997) corrected 
according to the prescription of North (1998), or, when no 
Geneva photometry was available, in the Stromgren system, 
applying the calibration of Moon & Dworetsky (1985), re- 
vised by Napiwotski et al. (1993). The luminosities have 
been obtained by taking into account the bolometric correc- 
tions measured by Lanz (1984). In case of binary systems, a 
correction for the duplicity has been applied to their magni- 
tudes (i.e., 0.75 mag for SB2 and 0.3 mag for SB 1 systems). 
The Lutz-Kelker correction (Lutz & Kelker 1973) has not 
been used, as it still remains a matter of debate (e.g., Arenou 
& Luri 2002). The stellar masses are obtained by interpola- 
tion in the evolutionary tracks of Schaller et al. (1992) for a 
solar metallicity Z = 0.018, as explained in North, Jaschek 
& Egret (1997). The radii have been computed from the 
luminosity and effective temperature. Finally, the surface 



© 2007 WILEY- VCH Verlag GmbH & Co. KGaA, Weinheim 



www.an-journal.org 



Astron. Nachr. / AN (2007) 



3 



Table 1 Fundamental parameters of stars with M > 3 M© . 
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gravity was obtained from mass and radius through its fun- 
damental definition. The errors on T e g and L determination 
were linearised and carefully propagated to errors on mass, 
radius, and log g, taking into account the slope of the evo- 
lutionary tracks. More details on the determination of the 
fundamental parameters can be found in our previous paper 
(Hubrig, North & Mathys 2000). 

The basic data for both sub-samples are presented in Ta- 
bles Q] and |2] The columns are, in order, the HD number of 
the star, the absolute visual magnitude, the mass (in solar 
masses), the logarithm of the effective temperature, the log- 
arithm of the luminosity (relative to the Sun), the logarithm 
of the stellar gravity, the radius (in solar radii), the distance 
d, the relative uncertainty of the parallax, and the fraction 
/ of the main-sequence life completed by the star, inter- 
polated from theoretical evolutionary tracks (Schaller et al. 
1992). For a number of stars located below the ZAMS the 
masses have been calculated by assuming their position on 
the ZAMS. Since ages are ill-defined for stars close to the 
ZAMS, in many applications we used gravity as a proxy 
for the relative age which has also the advantage of being a 



more directly measured quantity. For stars located close to 
the terminal-age main sequence (TAMS) we estimated an 
ambiguity of the order of 5% on the interpolated mass as 
there exists the possibility that these stars are making a loop 
back to the blue. 

3 Analysis and results 

3.1 The distribution of the magnetic Ap and Bp stars 
in the H-R diagram 

The position of the studied magnetic Ap and Bp stars in 
the H-R diagram is shown in Fig. Q] From this figure it is 
quite clear that the magnetic stars with M < 3 M Q show 
a concentration towards the centre of the main-sequence 
band and are rarely found close to either ZAMS or TAMS, 
in full agreement with our previous results (Hubrig, North 
& Mathys 2000, Hubrig, North & Szeifert 2005). By con- 
trast, the massive Bp stars are much younger, with a distri- 
bution crowded towards the ZAMS . We note that a number 
of stars in the mass domain between 3 and 4 M Q are dis- 
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Table 2 Fundamental parameters of stars with M < 3 M Q . 
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Fig. 1 H-R digram for Ap and Bp stars with accurate Hip- 
parcos parallaxes and well studied periodic variations of 
their mean longitudinal magnetic fields. The vertical and 
horizontal bars indicate the accuracy of the determination 
of effective temperature and gravity. 




Fig. 2 Distribution of relative ages for stars with masses 
M > 3 M Q (dotted line) and M < 3 M Q (solid line). 



tributed well across the main-sequence band. The distribu- 
tion of relative ages for each group is presented in Fig. [2] 
and shows a broad maximum indicating that the majority of 
magnetic stars with M < 3 Mq have completed 40 to 90% 
of their main-sequence life. The situation is completely dif- 
ferent for stars with M > 3 M : the distribution peaks at the 
ZAMS age with two secondary lower peaks at the relative 
ages around 60% and 80%. However, given the small num- 
ber of stars in this sample, the conclusion about secondary 
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Fig. 3 Cumulative distribution of log^-values for mag- 
netic stars with masses M > 3M Q (dashed line) and 
M < 3 M (solid line). 



peaks can be considered as only marginally significant and 
should be confirmed by larger statistics. 

Out of our sample of 57 stars with M < 3M Q , only 
three Ap stars, HD 12447, HD 32633, and HD 143473, have 
relative ages between 20 and 26%. In the following we show 
the difficulty of determining their fundamental parameters. 
HD 12447 with a relative age of 25% is a visual double star 
with Am = 1. The (B-V)T yc ho for this star was transformed 
to the (B-V) Johnson through the relation (B-V) JolmS on = 0.850| 
(B-V) T ycho given in ESA (1997), then to (B-V) G onova through| 
the relation given by Meylan & Hauck (1981). The sec- 
ondary has an index typical of an A5V star, the average 
colours of which (Hauck 1994) were subtracted from the 
colours of the binary to obtain the colours of the primary. 
Then T c g was obtained from X and Y and corrected ac- 
cording to the scheme proposed by North (1998). The T c g 
determination for HD 32633 with a relative age of 26% is 
very uncertain, ranging from our photometric estimate of 
12,819Kto 1 0,250 K using the (B2-G) index (Sokolov 1999).| 
We assumed T e g = 10,500 K, which corresponds to the spec-| 
tral type B9 given by Borra & Landstreet (1980). The red- 
dening Ay = 0.16 was adopted following Lucke's maps, 
instead of the initial Ay = 0.411 obtained from the X,Y 
Geneva parameters. The star HD 143473 with a relative age 
of 20% has already been discussed in our previous paper 
(Hubrig, North & Mathys 2000). The reddening obtained 
from the Geneva X, Y parameters appears much too high 
(Ay = 0.95) in view of the modest distance, and the effec- 
tive temperature for this star was therefore determined from 
spectroscopy rather from than from photometry. This star is 
also an extreme photometric Ap star with a photometric pe- 
culiarity parameter A(V1 - G) = 0.048 (Hauck & North 
1981). 

In Fig. [3] we show the cumulative distribution of \ogg 
for both sub-samples. A Kolmogorov-Smirnov test implies 
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Table 3 Model parameters for stars with M > 3 M Q . 



HD 


(Bi) 


frot 


v sini 


r 


i 


P 


5737 


324 


21.65 


20 


-0.583 


88.51 


5.64 


12767 


242 


1.89 


45 


-0.743 


33.32 


84.46 


18296 


440 


2.88 


25 


-0.618 


21.57 


84.67 


19832 


315 


0.73 


110 


-0.710 


47.00 


79.69 


21699 


828 


2.48 


35 


-0.803 


29.91 


86.40 


22470 


733 


0.68 


60 


-0.819 


14.59 


88.52 


24155 


1034 


2.53 


35 


-0.280 


63.82 


41.15 


25823 


668 


4.66 


25 


-0.214 


51.30 


51.05 


28843 


344 


1.37 


100 


-0.670 


89.92 


0.39 


34452 


743 


2.47 


46 


-0.371 


86.83 


6.88 


40312 


223 


3.62 


49 


-0.619 


44.74 


76.87 


49333 


618 


2.18 


65 


-0.780 


85.34 


33.40 


73340 


1644 


2.67 


- 


0.307 


- 


- 


79158 


672 


3,83 


60 


-0.910 


66.13 


83.92 


92664 


803 


1.67 


66 


0.163 


61.87 


21.04 


103192 


204 


2.34 


72 


0.301 


56.35 


19.68 


112413 


1349 


5.47 


10 


-0.670 


24.27 


84.91 


122532 


665 


3.68 


35 


-0.919 


56.57 


86.34 


124224 


572 


0.52 


119 


-0.510 


34.67 


77.35 


125823 


469 


8.82 


15 


-0.766 


48.33 


81.53 


133652 


1116 


2.30 


69 


-0.703 


84.36 


29.52 


133880 


2415 


0.88 


- 


-0.848 


- 


- 


137509 


1021 


4.49 


28 


-0.689 


50.72 


77.31 


142301 


2104 


1.46 


58 


-0.430 


45.16 


68.16 


142990 


1304 


0.98 


125 


-0.438 


54.60 


61.19 


144334 


783 


1.49 


44 


-0.194 


29.27 


69.28 


147010 


4032 


3.92 


22 


0.460 


71.24 


7.16 


151965 


2603 


1.61 


105 


0.216 


88.85 


0.74 


168733 


815 


14.78 


12 


0.306 


58.29 


18.17 


173650 


326 


9.98 


15 


-0.466 


45.19 


69.86 


175362 


3570 


3.67 


15 


-0.614 


25.80 


83.41 


196178 


973 


1.92 


50 


0.089 


51.33 


33.80 


223640 


643 


3.74 


30 


-0.064 


65.34 


27.56 



that the distribution of log g values for more massive stars 
differs from the distribution of stars with M < 3 M Q at a 
significance level of 99.8%. 

3.2 The evolution of the magnetic field strength across 
the main sequence 

Because of the strong dependence of the longitudinal mag- 
netic field on the rotational aspect, its usefulness to char- 
acterise actual magnetic field strength distributions is lim- 
ited, but this is overcome by repeated observations to sam- 
ple various rotation phases, hence various aspects of the 
magnetic field. As an estimate of the magnetic field strength 
we used the rms longitudinal magnetic field strength com- 
piled by Bychkov, Bychkova & Madej (2004). It has been 
computed from all measurements following the formula of 
Borra, Landstreet & Thompson (1983): 

/, " \ 1/2 

<^>=(^X>>* 2 J • (2) 



Table 4 Model parameters for stars with M < 3 M Q . 



TTT\ 

HL) 
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-0.357 


2.42 
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58.72 


1 1503 
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1.61 


13 


-0.479 
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87.71 
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62.00 
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56 


-0.810 
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14437 
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-0.934 


43.79 


88.13 


65339 


3113 


8.03 


12.5 
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0.660 


70.25 


4.21 


81009 


1431 


33.96 
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83368 
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r\r\r\ a a 

90044 


739 


4.38 


15 
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-0.096 


36.82 


58.31 


119419 


1770 


2.60 


34.8 
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-0.362 


30.69 


74.46 


152107 


1487 


3.87 


35 
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Fig. 4 The rms longitudinal magnetic field as a function 
of the completed fraction of main-sequence life. In this and 
the following figures, filled circles indicate stars with mass 
M > 3 Mo, while open circles indicate stars with mass 
M < 3M Q . 
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Fig. 6 The rms longitudinal magnetic field as a function 
of mass. 
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Fig. 7 Distribution of log P for stars with masses M > 
3 M (dotted line) and M < 3 M Q (solid line). 



Fig. 5 The rms longitudinal magnetic field as a function 
of effective temperature. 



The (B{) values are presented for each sub-sample sepa- 
rately in Tables[3]and|4]in the second columns and have been 
preferentially selected from magnetic field measurements 
carried out using Balmer hydrogen lines. As hydrogen is ex- 
pected to be homogeneously distributed over the stellar sur- 
face, the longitudinal magnetic field measurements should 
sample the magnetic field fairly uniformly over the observed 
hemisphere. In Fig.|4]the rms longitudinal magnetic field is 
plotted against the age of the stars expressed as a fraction of 



their total main-sequence life. For the sub-sample of mag- 
netic stars with M > 3 Mq the strongest magnetic fields 
are found in younger stars in terms of the elapsed fraction 
of their main-sequence life. The fact that the strongest mag- 
netic fields are only observed close to the ZAMS can be 
interpreted as a magnetic field decay in stars at advanced 
ages. A similar decay of the magnetic field in stars with 
M > 3 M Q had been suggested more than twenty years ago 
by North & Cramer (1984) from the study of logg values 
and surface magnetic fields estimated from Geneva photom- 
etry (see their Fig. 6a). However, even though their result 
was based on hundreds of stars, it had only a statistical and 
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Fig. 8 Rotation period as a function of mass. 




Fig. 9 Rotation period as a function of surface gravity for 
stars with M > 3M Q . The solid lines represent the the- 
oretical evolution of the period for a 4 M Q star and initial 
periods of 0.5 and 4 days, resulting from conservation of 
angular momentum alone. 



Fig. 10 Rotation period as a function of surface gravity 
for stars with M < 3 M Q . The solid lines represent the the- 
oretical evolution of the period for a 2 M Q star and initial 
periods of 0.5 and 4 days, resulting from conservation of 
angular momentum alone. 
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Fig. 11 Rotation period versus elapsed time on the main 
sequence. 



qualitative value, because the correlation they used between 
a photometric parameter and the surface magnetic field is 
not a one-to-one relation. On the other hand, the stellar radii 
of stars increase by approximately a factor of two during the 
stellar life time across the main sequence. Therefore, rather 
than a real decay of the magnetic field, our observations may 
simply testify for conservation of the magnetic surface flux: 
in such a case, the magnetic field strength is expected to 
have decreased by a factor of 4 by the time the star reaches 
the TAMS. 



The magnetic field for the majority of the stars with 
M < 3M0 becomes observable only after they have al- 
ready completed a significant part of their life on the main 
sequence. We note that the magnetic field strength obvi- 
ously decreases after the low mass stars have finished 80% 
of their main-sequence life. This would explain the rareness 
of detection of magnetic stars close to the TAMS. 

In Figs. [5] and [6] we present the rms longitudinal mag- 
netic field as a function of effective temperature and mass, 
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Fig. 12 The rms longitudinal magnetic field as a function 
of period. 



Fig. 14 Magnetic flux against elapsed time on the main 
sequence. 
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Fig. 13 Magnetic flux against rotation period. 

respectively. The distribution of stars in these diagrams re- 
veals that the strongest magnetic fields tend to be found in 
stars with an effective temperature between 1 0,000 K and 
15,000 K, and in the mass domain between 2.5 and 4 M Q . 

3.3 Rotation 

As mentioned above, our sample includes exclusively Ap 
and Bp stars with well studied periodic magnetic field vari- 
ations. The periods of variation are available in the cata- 
log of Bychkov, Bychkova & Madej (2005) and present a 
compilation from the literature and in some cases their own 
determinations. The rotation periods P for each sub-sample 
are listed separately in Tables [3] and [4] in the third column. 



Our previous study of the evolution of magnetic fields in 
Ap stars with well-known surface magnetic fields (Hubrig, 
North & Mathys 2000) suffered from a selection effect in 
the sense that the sample contained a high fraction (about 
2/3) of stars with rotation periods longer than 10 days, while 
the majority of the periods of magnetic stars fall between 2 
and 4 days. For the present study the distribution of rotation 
periods is representative of that of all Ap and Bp stars. In 
Fig. [7]we show the distribution of the periods for all stars in 
our sample. The majority of the studied stars have periods 
shorter than 4 days. Interestingly, the longest rotation peri- 
ods are found in stars in the mass domain between 1.8 and 
3 Mo (Fig. [8]l and with \ogg values ranging from 3.80 to 
4.13 (Figs. l9l and [Toll. In Fig. [TJJwe show the rotation pe- 
riod as a function of the relative age of Ap and Bp stars. The 
longest periods are found only in stars which spent already 
more than 40% of their main sequence life. For Bp and Ap 
stars with M > 3M Q , both Figs. |9"1 and [TT1 display a slight 
increase of rotation period with age which is consistent with 
the assumption of conservation of angular momentum dur- 
ing their life on the main sequence, without any hint of a 
braking mechanism. A similar result had already been ob- 
tained by North (1985) from the study of smaller samples of 
Ap and Bp stars with spectroscopic and photometric estima- 
tions of log g-values. North (1998) confirmed this trend in 
much the same way as here, except that his sample was not 
selected on the basis of magnetic fields, but on the basis of 
a confirmed Bp or Ap classification, photometric (or other) 
period and Hipparcos parallaxes. The observed strength of 
magnetic fields in stars with various rotation periods is pre- 
sented in Fig. [12] This figure clearly shows that the strongest 
magnetic fields are found only in fast-rotating stars. All stars 
with periods exceeding 10 days have (Bi) values less that 
2400 G. 
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In Figs. Qj] and [14] we present the magnetic flux '3/ (= 
12ttR 2 (Bi), assuming B s 3(5/)), plotted against the 
stellar rotation period and against the elapsed time on the 
main sequence. Our Fig. Qj] looks quite similar to Fig. 10 
presented in our previous study of magnetic stars with well- 
determined surface magnetic fields (Hubrig, North & Mathys| 
2000). Although the stars with M > 3 M show little signs 
of correlation between magnetic flux and rotational period 
— although one does have the visual impression of a slight 
positive correlation — the stars with M < 3 M Q do show 
some trend. The flux for the latter stars seems to take ran- 
dom values for logP < 1.2, while there is a slight de- 
creasing trend for logP > 1.5. It is also remarkable that 
there is a lack of points in the quadrant logP > 1.2 and 
log ^ < 26.7. The significance of this gap is attested by a 
2x2 contingency table, using the limits just mentioned: the 
two-sided Fischer exact probability value is p — 0.0251, so 
we are comfortably above the 95% significance level. On 
the other hand, an observational bias is not completely ex- 
cluded, because of the very long periods involved: it is pos- 
sible that the observers were motivated to follow only those 
stars with the largest magnetic fields for such long periods 
of time (remember that we selected here only stars with a 
decent magnetic curve). Thus the reality of the gap remains 
to be proven without any doubt. The plot of log \& against 
the elapsed fraction of the main sequence life shows signif- 
icant dispersion and no sign of evolution of the magnetic 
flux. 



3.4 Cluster members with accurate Hipparcos 
parallaxes 

Another way to get an insight into the evolution of stellar 
magnetic fields is to study magnetic stars in open clusters 
or associations at different ages. As far as the membership 
of Ap stars in distant open clusters or associations is con- 
cerned, we should keep in mind that such studies are mostly 
based upon photometry and upon radial velocity determi- 
nations. However, criteria for assessing cluster membership 
based on photometry cannot be applied to peculiar stars 
straight away, since strong backwarming effects lead to an 
anomalous energy distribution, thus affecting the position 
of the stars in colour-magnitude diagrams. 

In our sample, seven Ap/Bp stars are known members 
of nearby open clusters of different ages and have very ac- 
curate Hipparcos parallaxes (see Tabled- However, among 
them, only two stars, HD 92664 and HD 108945 have well- 
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Fig. 15 H-R digram for cluster stars. The vertical and hor- 
izontal bars indicate the accuracy of the determination of 
effective temperature and gravity. 



strong deficit of He. In summary, we are left with four clus- 
ter members which seem to represent typical magnetic stars. 

For most cluster stars, we preferred to use the distance 
given by main sequence fitting over the Hipparcos paral- 
lax. The error is generally assumed to be ±0.10 magni- 
tudes on the distance modulus. For the Pleiades, of which 
HD 23408 and HD 23950 are members, the adopted dis- 
tance was 132 ± 4pc, according to the geometric estimate 
by Zwahlen et al. (2004) for the Atlas binary, in agreement 
with the estimate of Munari et al. (2004) for the eclipsing 
binary HD 23642 (132 ± 2pc). For NGC 2451, of which 
HD 63401 is a member, we have adopted the distance mod- 
ulus of Maitzen & Catalano (1986): (m-M) = 6.7 and no 
reddening. For IC 2391, of which HD 74196 is a member, 
and for IC 2602, the true distance moduli were adopted from 
North (1993). For HD 108945 (in Coma Berenices), we just 
used the Hipparcos parallax, because of a possible depth ef- 
fect. 

The fundamental parameters for the cluster stars are pre- 
sented in Table [5] and their position in the H-R diagram 
is shown in Fig. [15] Both HD 92385 and HD 108945 with 
M < 3M spent already 30 and 80% of their main se- 
quence life, respectively, whereas HD 63401 and HD 92664 
with M > 3 M Q are very young and with a completed frac- 



defined magnetic rotational phase curves (Bychkov, Bychkovajion of their main sequence life of respectively 7 and 16%. 
& Madej 2005). The magnetic field measurements with FORS l|hese results are in full agreement with conclusions made 



at the VLT allowed to detect a magnetic field of the or- 
der of — 600 G in the spectra of two other stars, HD 63401 
and HD 92385 (Hubrig et al. 2006). On the other hand, we 
failed to measure magnetic fields in the remaining three 
stars, HD 23408, HD 23950 and HD 74196. It is quite pos- 
sible that they are not typical magnetic Ap/Bp stars, as they 
show in their spectra an overabundance of Mn and Hg and 



in the previous section in relation to the age of low and high 
mass stars. 

Only a few double-lined spectroscopic binary systems 
(which can be considered as very small clusters) containing 
a magnetic Ap star are currently known, and in all studied 
systems the Ap components with M < 3 M Q have already 
completed a significant fraction of their main-sequence life 
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Table 5 Fundamental parameters for cluster stars with accurate Hipparcos parallaxes. 



HD 


M v 


M/M Q 


log T cff 


log L/Lq 


log 9 


i?/R© 


d[pc] 


a(n)/iv 


/ 


23408 


-1.93 


4.463±0.104 


4.079±0.013 


2.916±0.053 


3.44±0.07 


6.66±0.57 


132. 


0.040 


1.11 


23950 


0.07 


3.401±0.103 


4.112±0.013 


2.167±0.053 


4.20±0.07 


2.42±0.21 


132. 


0.040 


0.26 


63401 


-0.37 


3.480±0.170 


4.129±0.013 


2.168±0.103 


4.28±0.10 


2.24±0.30 


219. 


0.109 


0.07 


74196 


-0.40 


3.604±0.114 


4.108±0.013 


2.352±0.056 


4.03±0.08 


3.05±0.27 


153. 


0.046 


0.57 


92385 


0.75 


2.730±0.106 


4.043±0.013 


1.804±0.077 


4.19±0.09 


2.19±0.23 


151. 


0.076 


0.30 


92664 


-0.32 


3.863±0.147 


4.154±0.013 


2.357±0.075 


4.24±0.09 


2.47±0.26 


143. 


0.073 


0.16 


108945 


0.57 


2.431±0.085 


3.944±0.015 


1.725±0.079 


3.82±0.09 


3.16±0.36 


95. 


0.079 


0.80 



(e.g., Wade et al. 1996, North et al. 1998, Carrier et al. 
2002). 

3.5 The magnetic field geometry 

There are two different ways to infer the geometry of the 
magnetic fields. Obviously, the observed magnetic curves 
can be modeled by assuming that the magnetic geometry 
can be described by the superposition of a dipole and a 
quadrupole magnetic field or higher multipoles with arbi- 
trary orientation. However, we should keep in mind that 
generally different magnetic field distributions and different 
inclination angles i between the rotation axis and the line 
of sight and inclination angles (3 between the magnetic axis 
and the rotation axis can represent the observed magnetic 
curves equally well. On the other hand, the vast majority of 
the magnetic CP stars exhibits a smooth, single-wave varia- 
tion of the longitudinal magnetic field during the stellar ro- 
tation cycle. These observations are considered as evidence 
for a dominant dipolar contribution to the magnetic field 
topology. A first statistical investigation of the orientation 
of magnetic axes in periodic magnetic variables based on 
the assumption that the studied stars are oblique dipole ro- 
tators has been carried out by Preston (1967) with the result 
that their magnetic and rotation axes are nearly orthogonal. 
He defined 

(B{) (min) cos (3 cos i — sin (3 sin i 



(Bi) (max) cos (3 cos i -t- sin [3 sin i ' 



so that 



f3 = arctan 



1 



cot i 



(3) 



(4) 



Observations of spectral line widths and profiles allow 
to determine the projected rotation velocity v sin i, where 
v is the equatorial rotation velocity and i is the inclination 
angle of the rotation axis. The equatorial rotation velocity 
is given by v cq = 50.6R/P, where R is the stellar radius 
in solar units and P the period in days. From the measured 
v sin i values and t> oq determined using the known rotation 
period and the radius computed from the luminosity and ef- 
fective temperature, the inclinations of the stellar rotation 
axes can easily be calculated. Knowing i and r, the de- 
termination of the obliquity angle (3 for the orientation of 
the magnetic axis with respect to the rotation axis is rather 
straightforward. 



For our study, the quantities r for all sample stars, ex- 
cept for HD 81009 have been taken from the catalogue of 
stellar magnetic rotational phase curves compiled by By- 
chkov, Bychkova & Madej (2005). For HD 81009 we used 
r = 0.64, determined by Landstreet & Mathys (2000). The 
v sin i values have been measured on our spectra or taken 
from published spectroscopic studies (Abt 2001, Abt, Lev- 
ato & Grosso 2002, Bagnulo et al. 2002, Balona & Laney 
2003, Bohlender, Landstreet & Thompson 1993, Hensberge 
et al. 1979, Landstreet 1970, Landstreet & Mathys 2000, 
Levato et al. 1996, Nielsen & Wahlgren 2002, Royer et al. 
2002, Wade et al. 2000). We could not find in the litera- 
ture values of vsini for three stars: HD54118, HD 73340 
and HD 133880. The inclination angles i between the rota- 
tion axis and the line of sight determined from v sin i and 
v cq have been used to calculate the obliquities following 
equation H] 

In columns 4 to 7 of tables [3] and |4j we present for 
each star in our sub-samples respectively the values v sin i, 
the parameter r from the catalog of Bychkov, Bychkova 
& Madej 2005 and the inclinations i and (3. For magnetic 
Ap stars with rotation periods of the order of one month 
or longer, the computed v cq are very small and the inclina- 
tion angle i between the rotation axis and the line of sight 
cannot be determined. A couple of years ago Landstreet & 
Mathys (2000) carried out modeling of observed magnetic 
field moments for a small sample of slowly rotating mag- 
netic stars. Their result was rather unexpected, implying that 
the magnetic stars with periods P>25 days have small val- 
ues of obliquity (3 of the model magnetic axis to the rotation 
axis, of the order of 20°. A few stars from their study, which 
have accurate Hipparcos parallaxes and available Stromgren 
and/or Geneva photometry have been included in our sam- 
ple. In fact, we simply employed the magnetic field param- 
eters determined by these authors in our study to increase 
the statistical significance of our results. 

In the upper panel of Fig. [16] we present the distribu- 
tion of r values for stars of M > 3 M and of M < 3 M Q 
shown by a dotted and a solid line, respectively. The distri- 
butions look rather different with an approximately bimodal 
distribution for lower mass stars with r peaked between 
—0.75 and —1 (highest peak) and between 0.25 and 0.5. 
The highest peak in the distribution of r values for stars with 
M > 3 Mq is shifted to higher values between —0.75 and 
—0.5. Both distributions nearly resemble the r value dis- 
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Fig. 17 r-value versus rotation period. 
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Fig. 16 Upper panel: Distribution of r-value for stars with 
masses M > 3 M (dotted line) and M < 3 M (solid line). 
Middle panel: r value versus mass. Lower panel: r-value 
versus effective temperature. 
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Fig. 18 The r-value versus elapsed time on the main se- 
quence. 



tribution presented in older studies of Preston (1967, 1971) 
and Landstreet (1970), with an apparent deficit of stars with 
intermediate values of (3. According to probability distribu- 
tions of r calculated by Landstreet (1970), we conclude that 
the majority of stars in our sample must have a (3 larger 60°. 
The lack of stars with r-values near zero can be explained 
either by an observational effect or by a lack of small (3- 
values. Further, we note that no star with M > 3M Q has 
r larger than 0.5, indicating the difficulty of detecting pe- 
riodicity in stars with very small magnetic field variations. 
The examination of the middle panel of Fig. [16] with r val- 
ues plotted against stellar mass reveals no obvious depen- 
dence for lower mass stars whereas higher mass stars with 
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Fig. 19 Distribution of obliquity angles (3 for stars with 
masses M > 3 M Q (dotted line) and M < 3 M Q (solid line). 
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Fig. 21 /3-value for high (upper panel) and low (lower 
panel) mass stars versus elapsed time on the main sequence. 



Fig. 20 (3 versus rotation period. 

M > 4M0 show a slight trend to have r < —0.5. This 
behaviour can be interpreted by the prevalence of larger 
obliquities (3 in more massive stars. From the inspection of 
the dependence between the r-values and effective temper- 
atures (lower panel in Fig.fToT) the same is possibly true for 
the hottest stars in our sample, i.e. larger obliquities (3 are 
found in hotter stars. 

Quite interesting results follow from the study of the dis- 
tribution of r-values for stars with different rotation periods 
and the distribution of r-values over the age, i.e. elapsed 
time on the main sequence presented in Figs. [T7] and [18] 
While we observe a random distribution of obliquity angles 
(3 in faster rotating stars, it seems that long period magnetic 



variables of stars with M < 3 M Q tend to have values of 
(3 close to 0°. In the distribution of r-values over the age 
we observe that the obliquity angle distribution as inferred 
from the distribution of r-values appears random at the time 
magnetic stars become observable on the H-R diagram. Af- 
ter quite a short time spent by a star with mass M < 3 M Q 
on the main sequence the obliquity angle (3 tends to reach 
values close to either to 90° or 0°. In other words, the mag- 
netic axis becomes either perpendicular to the rotation axis 
or aligned with it with advanced age. A similar result re- 
vealing the bimodal distribution based on a smaller sample 
of Ap and Bp stars had been presented by North (1985), but 
on the basis of less reliable photometric log g estimates. For 
more massive stars with M > 3 M Q the obliquity angle (3 
tends to reach values close to 90°. 
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Fig. 22 /3-value versus mass (upper panel) and /3-value 
versus effective temperature (lower panel). 



The distribution of obliquity angles j3 for stars with M > 
3M and M < 3M is presented in Fig. [T9j In agree- 
ment with the discussion on the distribution of r-values, 
the excess of stars with large /3-values is detected in both 
higher and lower mass stars. Recently, Bychkov, Bychkova 
& Madej (2004) presented a study of parameters of mag- 
netic phase curves based on their own compilation in By- 
chkov, Bychkova & Madej (2005) of the data for 1 34 Ap and 
Bp stars. Surprisingly, they found an excess of stars with an- 
gles (3 close to 0° . We do not know the origin of this discrep- 
ancy between their and our study, but we only emphasise 
here that our sample consists of carefully selected 90 Ap and 
Bp stars with a sufficient number of magnetic field measure- 
ments and well-defined periods and magnetic phase curves. 
From the inspection of plots displaying the distribution of (3 
in fast and slowly rotating stars (Fig. [20b. it is obvious that 



large j3 (>80°) outweigh the distributions of faster rotating 
stars whereas the angle (3 is smaller than 20° in slower ro- 
tating stars. We note that all nine points with small values of 
/3 in slowly rotating stars are from the study of Landstreet 
& Mathys (2000). As for stars with M > 3M , the two 
stars with the longest periods (HD5737 and HD 168733) 
have magnetic and rotation axes aligned to within 20°, sim- 
ilar to the behaviour of slowly rotating stars of lower mass. 
Clearly, more magnetic field studies of slowly rotating stars 
with higher mass are needed to confirm this behaviour. 

The evolution of the obliquity angle j3 seems to be quite 
different for low and high mass stars (Fig. I2TI 1. We find a 
strong hint for an increase of (3 with elapsed time on the 
main sequence for stars with M > 3 M . However, no sim- 
ilar trend is found for lower mass stars, although the pre- 
dominance of high values of (3 at advanced ages is notable. 
Further, we do not find any trend for /3-values either with 
mass or with T e g (Fig.l22l). 

As is mentioned by Bychkov, Bychkova & Madej (2005),| 
a number of stars in their catalog exhibit marked anhar- 
monicity, which indicates significant departures from dipole 
geometry. To reproduce the shapes of the observed longitu- 
dinal magnetic field variation, the magnetic field distribu- 
tion has to be regarded as a sum of dipole and quadrupole 
components. Out of the 33 stars studied with M > 3M , 
21% have magnetic phase curves fitted by a double wave, 
whereas the percentage of stars with M < 3 M for which 
the magnetic field variations had to be fitted by a double 
wave is 12%. This statistics could indicate the possibility 
that the magnetic topology in stars with M > 3 M is fre- 
quently more complex than that in less massive stars. The 
double-wave longitudinal field curves in stars with M < 

3 M are mostly found at the age of 40-50% of their main 
sequence life, while a number of stars with M > 3 M ex- 
hibit a non-dipolar field already at a very young age close 
to the ZAMS. However, further systematic studies of mean 
longitudinal magnetic fields should be conducted to better 
assert the evolution of the magnetic field distributions dur- 
ing the main-sequence life. Moreover, the observations of 
other magnetic field moments such as the mean quadratic 
field and mean crossover field, along with studies of Zee- 
man structure in the individual line profiles in circular and 
linear polarisation should be involved in future modeling to 
obtain a most accurate determination of the magnetic field 
topologies. 

4 Discussion 

We believe that the work presented here has important im- 
plications for the understanding of the origin of the mag- 
netism detected in A and B stars. The new results based on 
the statistical properties of 90 Ap and Bp stars can be sum- 
marised as follows: 

- The present study of the evolutionary state of upper main| 
sequence magnetic stars indicates a notable difference 
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between the distributions of high mass and low mass 
magnetic stars in the H-R diagram. In contrast to mag- 
netic stars of mass M < 3 M© which are mostly found 
around the centre of the main-sequence band, the stars 
with masses M > 3 M Q seem to be concentrated closer 
to the ZAMS, and the stronger magnetic fields tend to be 
found in hotter, younger (in terms of the elapsed fraction 
of main-sequence life) and more massive stars. 

- The longest rotation periods are found in stars in the 
mass domain between 1.8 and 3 M Q and with log g val- 
ues ranging from 3.80 to 4.13. The longest periods are 
found only in stars which spent already more than 40% 
of their main sequence life. 

- No evidence is found for any loss of angular momentum 
during the main-sequence life of stars with M > 3M Q . 

- No sign of evolution of the magnetic flux over the stel- 
lar life time on the main sequence has been found, in 
agreement with the assumption of magnetic flux conser- 
vation. 

- Our study of a few known members of nearby open clus- 
ters of different ages and with accurate Hipparcos par- 
allaxes confirms the conclusions based on the study of 
field stars with accurate Hipparcos parallaxes. 

- The excess of stars with large obliquities [3 is detected 
in both higher and lower mass stars. It is quite possible 
that the angle (3 becomes close to 0° in slower rotating 
stars of high mass too, analog to the behaviour of angles 
(3 in slowly rotating stars of M < 3 M Q , 

- The obliquity angle distribution as inferred from the dis- 
tribution of r-values appears random at the time mag- 
netic stars become observable on the H-R diagram. Af- 
ter quite a short time spent by a star on the main se- 
quence the obliquity angle (3 tends to reach values close 
to either 90° or 0° for M < 3 M©. However, for M > 
3 M the obliquity angle (i in increasing with the com- 
pleted fraction of the main-sequence life. It would be 
important to consider such a behaviour in the frame- 
work of dynamo or fossil magnetic theories, taking into 
account the influence of various mechanisms (e.g. possi- 
ble changes in internal rotation during the evolutionary 
stage). 

- While we find a strong hint for an increase of [3 with 
the elapsed time on the main sequence for stars with 
M > 3 M Q , the trend for lower mass stars is much less 
pronounced, though it seems to exist as well. 

The presented results have to be implemented in the 
long-lasting discussions about fossil or dynamo origin for 
the observed magnetic fields. It is quite possible that the ob- 
served magnetic fields in stars with M > 3M Q , showing 
the prevalence of stronger magnetic fields in younger stars, 
are in some sense fossil, being the remnants of magnetic 
fields originally present in the material from which these 
stars formed. Recent theoretical work indicates that dynamo 
action is possible, even likely, in a convective core, but the 
overlying stellar radiative envelope represents a significant 
impediment to the appearance of any of the generated mag- 



netic field at the surface (e.g., MacGregor & Cassinelli 2003,| 
MacDonald & Mullan 2004). To remedy the difficulties as- 
sociated with the transport of magnetic fields from the core 
to the surface there have been a number of attempts to iden- 
tify mechanisms that are capable of generating a magnetic 
field within the radiative envelope itself (e.g., Spruit 2002, 
Braithwaite 2004, Arlt, Hollerbach & Rudiger 2003). For 
stars with mass M < 3 M© the magnetic fields appear after 
they completed more than 20-30% of their main-sequence 
life. Because of this apparent difference in the evolutionary 
state between higher mass and older lower mass stars it is 
tempting to make a conclusion that magnetic fields in stars 
of lower mass are generated by dynamo action. However, 
the comparison of the observed magnetic field evolution and 
magnetic field geometries in both samples presented in this 
paper does not reveal any striking differences apart from 
the older age of lower mass magnetic stars on the main se- 
quence, the very fact that most extreme slow rotators are 
found among stars with M < 3M Q , and the difference in 
the evolutionary behaviour of obliquities between both sam- 
ples. 
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